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ABSTRACT 

We report the discovery of a powerful and variable off-nuclear X-ray source in the 
nearby spiral galaxy NGC 4945. Two ROSAT PSPC observations show the source to 
brighten in 0.5-2.0 keV flux by a factor of « 9 on a time-scale of 11 months or less. 
It is seen by ASCA about one month after the second PSPC pointing, and is seen to 
have dimmed by a factor of ^ 7 in a ROSAT HRI pointing about one year after the 
second PSPC pointing. Its maximum observed 0.8-2.5 keV luminosity is ^ 8 x 10'^* 
erg s~^, making it brighter than any known persistent X-ray binary in the Milky Way. 
Its total X-ray luminosity is probably larger than 1.2 x 10^^ erg s~^. The observed 
variability argues against a superbubble interpretation, and the off-nuclear position 
argues against a low-luminosity active galactic nucleus. The source is therefore prob- 
ably either an ultra-powerful X-ray binary or an ultra-powerful supernova remnant. 
Optical monitoring has not identified any supernovae in NGC 4945 during the time of 
the X-ray observations, and any supernova would have had to have been either very 
highly absorbed or intrinsically optically faint. 
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1 INTRODUCTION 

Recent X-ray studies of nearby spiral galaxies have revealed 
a wealth of activity (e.g. Fabbiano 1995 and references 
therein). Low- luminosity active galactic nuclei, starburst su- 
perbubbles of hot gas, and ultra-powerful X-ray binaries and 
supernova remnants are all currently subjects of widespread 
interest. Luminous and variable off- nuclear X-ray sources 
are of particular interest, as they often challenge our un- 
derstanding of accretion processes in stellar systems (e.g. 
CoUura et al. 1994; source A in IC 342, Okada et al. 1994). 

NGC 4945 is a nearby (3.7 Mpc: Mauersberger et al. 
1996), edge-on spiral galaxy in the southern hemisphere. It 
is the third brightest galaxy in the IRAS point source cata- 
logue, and most of its infrared emission arises in a compact 
central region (Rice et al. 1988; Brock et al. 1988). NGC 
4945 shows both starburst emission (e.g. Heckman, Armus 
& Miley 1990; Koornneef 1993) and a highly obscured and 
variable Seyfert core (Iwasawa et al. 1993). It has recently 
been discovered to be the brightest Seyfert 2 known at high 
X-ray energies (« 100 keV: Done, Madejski & Smith 1996). 
In this Letter we report on a powerful, off-nuclear X-ray 
source in NGC 4945 that shows large-amplitude variability 
on a time-scale of a year. 



2 OBSERVATIONS AND DATA ANALYSIS 

NGC 4945 was observed with the ROSAT Position Sensitive 
Proportional Counter (PSPC; Triimper 1983; Pfeffermann 
et al. 1987) starting on 1992 August 12 (RP600259: total 
raw exposure of 14208 s over 2.66 x lO'^ s) and 1993 July 10 
(RP600275: total raw exposure of 9032 s over 80760 s). These 
were the first observations of this galaxy with an imaging X- 
ray detector. Hereafter we shall refer to RP600259 as 'Point- 
ing 1' or 'PI' and RP600275 as 'Pointing 2' or 'P2'. The 
ROSAT observations were performed in the standard 'wob- 
ble' mode, and reduction and analysis of the PSPC data 
were performed with the Starlink ASTERIX X-ray data pro- 
cessing system. NGC 4945 was also observed with ASCA 
(Tanaka, Inoue & Holt 1994) starting on 1993 August 31 
(total good SIS exposure time of 35000 s over 68362 s), 
and the ROSAT High Resolution Imager (HRI) starting on 
1994 July 11 (RH600276: total raw exposure of 8608 s over 
6.27 X 10^ s). The ASCA data were reduced and analysed 
using the xselect X-ray data processing system, and the 
HRI data were reduced and analysed with asterix. 

In Figure 1 we show contours of the P2 PSPC hard 
image (we use PSPC channels 52-201 which correspond to 
«0. 5-2.0 keV) overlaid on the image from the UK Schmidt 
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NGC 4945 - optical image and PSPC hard band image 
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Figure 1. Contours of the P2 ROSAT PSPC X-ray image overlaid on the image from the UK Schmidt Southern Sky Survey J plate. 
Contours are at 8.8, 17.6, 35.2 and 70.6 per cent of the maximum pixel value. The location of the variable X-ray source, NGC 4945 X-1, 
is marked with an eight-pointed star. The nuclear position of Henkel at al. (1994) is marked with a solid white dot. Note that NGC 4945 
X-1 is significantly offset from the nucleus. 



Southern Sky Survey J plate [see section lib of Lasker et al. 
(1990) for more information on the optical image]. We also 
show the position of the nucleus (we use the nuclear co- 
ordinates of Henkel, Whiteoak & Mauersberger 1994). The 
ROSAT imago shows several point sources as well as diffuse 
X-ray emission (presumably associated with the starburst). 
We have used the ASTERIX point source searching program 
PSS (Allan 1995) to search the hard PSPC image for point 
sources that are spatially coincident with the optical extent 
of NGC 4945. Wo detect five point sources (at greater than 
5a) in both PI and P2. If we compare the count rates of these 
sources between PI and P2, four of them show little evidence 
for variability. However, the fifth stands out as showing clear 
evidence for having brightened between PI and P2 (see Fig- 



ure 2). We have plotted the position of this source in Figure 
1 and shall hereafter refer to it as NGC 4945 X-1. The PI 
X-ray centroid position of X-1 is a2ooo = 13*' 05™ 22.6'', 
(52000 = —49° 29' 12", and its P2 X-ray centroid position is 
02000 = 13^ 05™ 22.1% (52000 = -49°29'09". These positions 
each have absolute errors of about 15 arcsec and are entirely 
consistent with one another. X-1 is offset from the nucleus 
with extremely high statistical significance (the separation 
is Ri 1.3 arcmin). We have checked the PSPC image eistrom- 
ctry (for both PI and P2) using NGC 4945A, which lies 
close to NGC 4945 on the plane of the sky, and it appears 
to be good to within 15 arcsec. The PSS count rate of X-1 
during PI is (2.7 ±0.6) x 10~^ count s~\ and its PSS count 
rate during P2 is (24.7 ± 2.3) x 10~^ count s"^ (we quote 
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NGC 4945 X-1 light curve 
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Figure 2. Effective PSPC count rate versus time for NGC 4945 
X-1. The effective PSPC count rates from ASCA and ROSAT 
HRI data are computed as per the text. 



68 per cent confidence count rate errors). The hard ROSAT 
band count rate thus apjjcars to have increased by a factor 
of 9.1 ± 2.3 in the 11 months between PI and P2. We note 
that, while small-scale gradients in the diffuse X-ray emis- 
sion might lead to slightly larger errors than the formal PSS 
errors stated above, they cannot remove the large-amplitude 
variability that we detect. Inspection by eye of the PI and 
P2 images clearly shows that X-1 has brightened substan- 
tially. During P2, X-1 is the brightest soft X-ray point source 
observable in NGC 4945. 

The ASCA observation also detects a point X-ray source 
with a position that coincides with that of X-1. We have used 
other sources in NGC 4945 seen by both ASCA and ROSAT 
to improve the astrometry of the ASCA SIS images. The 
ASCA X-ray centroid position for the source is 02000 ~ 13*^ 
05™ 22.9% 52000 = -49°29'17", and the error circle for this 
position is about 20 arcsec. The source is detected with very 
high statistical significance (> 6cr). Precision source count 
rate determination is difficult, owing to the broader ASCA 
point spread function that leads to cross-source contamina- 
tion. In order roughly to estimate the count rate, we have 
compared the relative peak heights of X-1 and a source lo- 
cated at 02000 = 13'' 05"" 31.9% ^2000 = -49°27'33" in the 
0.6-2.0 keV ASCA image. The height of X-1 is larger than 
that of the comparison source by a factor of about 1.5. Pro- 
vided that the comparison source is not variable we estimate 
an effective PSPC count rate for X-1 of about 11 x 10~^ 
count s~^ by multiplying the PSPC count rate of the com- 
parison source by a factor of 1.5. Thus X-1 appears to have 
dimmed by a factor of about 2 in 52 days. We note that X-1 
is not the brightest point source in the 0.6-2.0 keV ASCA 
image of NGC 4945. 

We have also used PSS to search the ROSAT HRI image 
for point sources. In this searching we use only HRI channels 
4-16 to avoid contamination by ultraviolet light (cf. section 
3.6 of David et al. 1995). At the PSPC position of X-1 there 



is only a weak (l.Scr) X-ray source with an HRI count rate of 
(0.45 ±0.34) X 10"'^ count s"\ The HRI count rate of X-1 is 
thus a factor of ^ 28 times lower than its P2 count rate. If 
we use PIMMS (Mukai 1995) and the PSPC spectral models 
given below to take into account the different responses of 
the PSPC and the HRI, we find that the HRI count rate 
of X-1 should have only been lower than its P2 count rate 
by a factor of 2.5-4.0 if X-1 remained constant. The HRI- 
to-P2 count rate ratio of X-1 is much smaller than that for 
the other sources in NGC 4945. Thus X-1 appears to have 
dimmed by a factor of ^ 7 between 1993 July 10 and 1994 
July 11. 

Count rates from ROSAT should be averaged over an 
integer multiple of the 400-s wobble period when used to 
search for rapid variability of cosmic X-ray sources (cf. 
Brinkmann et al. 1994). We have searched for rapid vari- 
ability of X-1 during PI and P2 and do not detect any such 
variability with high statistical significance. However, oven 
during P2 the statistics are such that we would only be able 
to detect variability by factors ^ 2. The ASCA data do not 
show any significant X-1 variability (again variations by a 
factor of ^ 2 would be required for a detection), and the 
HRI data do not allow a serious variability search. 

Meaningful PSPC spectral fitting of X-1 is not possi- 
ble during PI, and only crude spectral fitting is possible 
during P2. Cross-source contamination, the low total num- 
ber of source counts, and the limited spectral resolution of 
the PSPC all hinder precision fitting. However, in order to 
probe the spectrum of X-1 to the greatest extent possible 
and derive luminosity constraints, we have extracted source 
counts from a circular region of radius 0.8 arcmin centred on 
X-1. In this region the emission from X-1 dominates the X- 
ray flux. While we recognize that we are losing some source 
counts as a result of our small extraction radius, any larger 
radius would introduce unacceptably large cross-source con- 
tamination. We do include the standard energy- dependent 
correction for source counts scattered outside our source re- 
gion. We take our background from a circular region of rar 
dius 1 arcmin centred on 02000 = 13*" 05™ 10. 8^^, ^2000 = 
— 49°31'26", as this region appears to have approximately 
the same level of background as the source region. We obtain 
186 ± 15 counts for our spectrum after background subtrac- 
tion. The Galactic neutral hydrogen column towards NGC 
4945 is (1.1 ±0.3) x 10^^ cm"^ (Heiles & Cleary 1979, page 
51), and we note that there may well be signiflcant intrinsic 
column in the disc of NGC 4945. 

After corrections for source region size, vignetting, wire 
scattering and detector dead time, we have binned our ex- 
tracted spectrum so that there arc at least 12 counts in 
each spectral bin (this allows the use of chi-squared fitting 
techniques). We fit this spectrum using the xspec spectral 
fitting package. We fit only the data above 0.85 keV, as the 
lower energy data have potentially significant background 
and cross-source contamination (by ignoring the data be- 
low 0.85 keV we lose only about 10 per cent of our source 
counts due to the large column). When we state fluxes be- 
low they are computed for the 0.85-2.5 keV band that we 
directly flt. All errors are for 90 per cent confldence, taking 
all free parameters to be of interest other than the abso- 
lute normalization. We flrst considered a power-law model 
with the only absorption being due to the Galactic col- 
umn. This model gives a flat power-law photon index of 
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r = 0.02 ± 0.7, xl = 0.67 for 12 degrees of freedom, an ab- 
sorbed flux of 4.8 X 10~^^ erg cm~^ s~^, and an unabsorbcd 
flux of 5.2 X 10~^^ erg cm"^ s~^. If we instead fix the pho- 
ton index to be F = 2.0 and allow the absorption column 
to be a free parameter, we obtain A^h = (7.0j^2 4) ^ ^0^^ 
cm~^, = 0.74 for 12 degrees of freedom, an absorbed 
flux of 4.0 X 10~^^ erg cm~^ s~^, and an unabsorbed flux 
of 8.0 X 10~^^ erg cm~^ s~^. Power-law models where we 
allow both the photon index and the absorption column to 
be free parameters arc very poorly constrained, but give un- 
absorbed fluxes similar to or larger than those stated above. 
We have also fltted an absorbed thermal bremsstrahlung 
model to our spectrum (sec Figure 3). Using this model we 
derive iVn = (e.Otg.a) x 10^^ cm"^ kT = 4.5 kcV [kT is 
constrained to be larger than 0.5 keV), ~ 0.80 for 11 de- 
grees of freedom, an absorbed flux of 4.1 x 10~^^ erg cm~^ 
s~^, and an unabsorbed flux of 7.1 x 10~^^ erg cm~^ s~^. 
If we fix the total column at the Galactic value, a thermal 
bremsstrahlung model leaves large systematic residuals and 
can be ruled out with greater than 90 per cent confidence. 

The ASCA cross-source contamination does not allow 
a detailed spectral analysis, but we have performed a rough 
X-ray colour analysis using a small aperture (a square with 
side length 1.2 arcmin) to minimize cross-source contamina- 
tion. The colour analysis used three energy bands: 0.7-1.3 
keV, 1.3-2.3 keV and 2.3-8.0 keV. X-1 has the smallest 2.3- 
8.0 keV to 1.3-2.3 keV count rate ratio (0.84 ± 0.15) as 
well as the smallest 0.7-1.3 keV to 1.3-2.3 keV count rate 
ratio (0.39 ± 0.07) of any source in NGC 4945. For a power- 
law model, the first ratio corresponds to a spectral slope of 
r = 1.7 or steeper, and the second ratio corresponds to an 
absorption column of 5 x lO^'^ cm~^ or more. This suggests 
that the Galactic column power-law fit to the PSPC data 
(see above) is probably not physically appropriate. A more 
detailed colour analysis of all the sources in NGC 4945 will 
be presented by Iwasawa et al. (in preparation). 



3 DISCUSSION 

It is unlikely that X-1 is a foreground or background X- 
ray source. Given the position of NGC 4945, the probability 
of finding a contaminating source with at least the P2 flux 
of X-1 within the optical extent of NGC 4945 is less than 
1 X 10"^ (see section V of Hertz & Grindlay 1984). We have 
searched the UK Schmidt plate image and the less saturated 
image in Sandage & Bedke (1994), and there axe no obvious 
foreground optical counterparts. The position of X-1 in the 
nomograph of Maccacaro et al. (1988) docs not agree well 
with that of a foreground star (see their Figure 1). There is 
no strong X-ray variability during P2 as might be expected 
from a foreground flare star and, in addition, ASCA saw X- 
1 to be bright 52 days after P2. Furthermore, we note that 
the excess column over the Galactic value suggested by our 
spectral analysis argues for the extragalactic nature of X-1. 

As all of our models above yield unabsorbed fluxes 
larger than 5 x lO^^'' erg cm~^ s"'^, the isotropic luminos- 
ity of X-1 during P2 appears to be larger than 8 x 10** erg 
s~^ (we use a conservative distance of 3.7 Mpc). Wc note 
that there is probably some additional flux outside the 0.85- 
2.5 keV band considered above, and reasonable estimates of 
the total X-ray flux of X-1 suggest isotropic luminosities 
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Figure 3. ROSAT PSPC spectrum of NGC 4945 X-1 during P2. 
An absorbed thermal bremsstrahlung model has been fitted to 
the data (see the text), and the data-to- model ratio is shown. 

of 1.2 X 10^^ erg s^^ or more. There is no obvious radio 
source that corresponds to X-1 in the 843 MHz radio map 
of Harnett & Reynolds (1985) or the CO maps of Dahlem 
et al. (1993), although these limits are not tight (e.g. the 843 
MHz radio flux is less than 77 mjy with a 43 x 57 arcsec'^ 
beam). The large amplitude variability of X-1 allows us to 
rule out a starburst superbubble model for this source, and 
its off'-nuclear position argues against a low luminosity active 
galactic nucleus. We are thus led to consider ultra-powerful 
X-ray binaries and supernova remnants. 

The minimum P2 isotropic luminosity for X-1 corre- 
sponds to the Eddington limit luminosity for a 6 solar 
mass compact object, and makes it more luminous than any 
known persistent X-ray binary in the Milky Way. If X-1 is 
a binary, then the variability and the mass derived from 
the Eddington limit are suggestive of a transient black hole 
binary, and we note that the peak observed isotropic lumi- 
nosity of X-1 is similar to that for the Galactic transient GS 
2023-F338 (V404 Cyg; Tanaka & Lewin 1995). However, we 
note that some rare neutron star binaries in the Magellanic 
Clouds have comparable isotropic luminosities to X-1 (e.g. 
A0535 - 668 in outburst; SMC X-1 in its bright state). Of 
course, the emission need not be isotropic and indeed could 
be beamed (see Reynolds et al. 1996 for a more detailed 
discussion of this possibility) . Our X-ray spectrum is consis- 
tent to within its (large) errors with that seen from luminous 
Galactic black hole binaries. 

Ultra-powerful supernova remnants arc described in 
Schlegel (1995) and references therein. To explain X-1 a rem- 
nant would have had to brighten by a factor of about 9 in one 
year and then dim by a factor of ^ 7 in the subsequent yeax. 
This could perhaps happen if we caught a remnant either 
when it was very young or when it happened to bo interact- 
ing with circumstellar material. By comparison with predic- 
tions for the X-ray emission from the ring around SN 1987a 
(e.g. Suzuki, Shigeyama & Nomoto 1992), we note that very 
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dense circumstellar material would be needed in the second 
scenario to explain the large luminosity and short cooling 
time. However, the required variability does not seem to 
be obviously inconsistent with our limited knowledge of the 
variability of other powerful supernova remnants (Schlogol 
1995). We have checked the updated version of the Asiago 
supernova catalogue and the van den Bergh supernova cat- 
alogue (Barbon, Cappellaro & Turatto 1989; van den Bergh 
1994), and no supernovae in NGC 4945 are listed. NGC 4945 
is regularly monitored for supernovae, and no supernovae 
brighter than 15.5 mag have been seen during the dates rel- 
evant to this paper (Rev. R. Evans, private communication; 
see chapter 5 of Murdin 1990). We note that the 'inclination 
effect' is not likely to make the discovery of supernovae im- 
possible in NGC 4945 (see the discussion in van den Bergh 
1994). In the case of only Galactic extinction, the apparent 
magnitude limit corresponds to an absolute magnitude limit 
of —12.9. If we take into account the intrinsic extinction sug- 
gested by our X-ray spectral fitting (see above), the absolute 
magnitude limit is reduced to about —15.6, although limit 
values as low as —20 are not impossible. However, we note 
that if the extinction were eis large as 7.1 magnitudes then 
even larger 0.85-2.5 keV luminosities of 1.7 x 10^^ erg s~^ 
or more would be implied. Thus if X-1 wore a supernova it 
probably, although not certainly, would have been seen at 
optical wavelengths. We have also checked the UK Schmidt 
plate database for optical plates that might be used to look 
for a supernova, but unfortunately there are no plates at 
appropriate dates. 

Additional ASCA, ROSAT and SAX observations of 
NGC 4945 are needed to monitor the variability of X-1. 
The detection of another outburst or rapid X-ray variability 
would further support the binary interpretation. 
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